NACA TN 8542 Qggg 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 3542 



9 




) 


ANALYSIS OF STRESSES IN THE PLASTIC RANGE AROUND 

A CIRCULAR HOLE IN A PLATE SUBJECTED 

TO UNIAXIAL TENSION 

By Bernard Budiansky and Robert J. Vidensek 

Langley Aeronautical Laboratory 
Langley Field, Va. 



Washington 
October 1955 



I tech library kafb, niw 


G 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECH UBRARY KAFB, NM 


TECHNICAL NOTE 35^2 


ANALYSIS OF STRESSES IN THE PLASTIC RANGE AROUND 
A CIRCULAR HOIE IN A PLATE SUBJECTED 
TO UNIAXIAL TENSION 

By Bernard Budlansky and RoLert J. Vidensek 




>t>3 


SUMMARY 


An approximate theoretical solution is presented for the stresses 
in the plastic range around a circular hole in an infinite sheet sub- 
jected to unlaxieJ. tension. The solution is based on the sin^jle defor- 
mation theoiy of plasticity and is found by appJjLcation of a variational 
principle in conjunction with the Rayleigh-Ritz procedure and the use of 
a high-speed con^iutlng machine (SEAC). Numerical results are obtained 
for four different materials, which are characterized by four distinct 
uniaxial stress-strain curves. The results for stress concentration 
factor in the plastic range are coinpared with those obtained from a for- 
mula due to Stowell. 


INTRODUCTION 


The stress distributions that occur around structural discontinuities 
such as holes and notches have been found theoretically for a wide variety 
of canes on the basis of the theory of elasticity. An important problem 
is the corresponding determination of such stress distributions for 
strain-hardening materials when the stresses exceed the elastic limit. 

A major obstacle to such an undertaking lies in the fact that basic 
stress-strain relations in the plastic range have not yet been definitely 
established for strain-hardening materials; in addition, even after the 
choice of a particular stress-strain relation is made, the concomitant 
nonlinear system of equations governing the stress distribution may 
generally be expected to defy exact analytical solution. 

The present paper considers the problem of finding the stresses in 
the plastic range around a clrciilar hole in an infinite plate subjected 
to uniaxial tension at infinity (fig. l). The plate material is assumed 
to obey the stress -strain relations of the sin 5 >le deformation theory of 
plasticity, and an approximate solution is effected by application of an 
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apprcJpriate variational principle in conjunction with the Rayleigh-Ritz 
procedure and the use of a high-speed con^juting machine, the Standards' 
Eastern Automatic Computer (SEAC) of the National Bureau of Standards. 
Numerical results are given for fovur different materials, each of which 
is characterized hy a particular uniaxial stress-strain curve. 

The results found for the stress concentration factor in the plastic 
range are conipared with those predicted hy a siirple formula suggested hy 
StoweU in 1950 (ref. l). 


SYMBOLS 


<^r-> '''^rG 


radial, circumferential, and shear stresses, 
respectively 


— o _ o o 


"^rG 


radial, circumferential, and shear stresses corre- 
sponding to elastic solution, respectively 

correction stresses (for example, Oj. = - aj.°) 

\l/2 


Cg effective stress, 

cTe° = (cr°f + (cT0°/ - a/ae° + 5(T^e°) 


A 2 - 2 

*^e “ \^r + "^0 


A = Uoc/ai 
K 

^r^’ ^0-> ^rG 


- 2A1/2 

'^r^G ^^rG / 


uniaxial stress at Infinity 


nominal yield stress in Ramiberg-OBgood eq^uatlon 
(see eq.. (3)) 


stress concentration factor 

radial, circumferential, and shear strains, 
respectively 


radial and circumferential displacements, 
respectively 


Ur, uq 
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t 

apq_ 


0 

P 

^ = l/p 
E 
Es 
G 

V 

n 

w 

F 

Emod 

0 

i 

j 

M, N 

q. 


radial and circumferen'tial displacemen'ts corresponding 
to elastic solution, respectively 

correction displacements (for exan 5 >le, Up = u^. - Up°) 
stress function for a^, Oq, n.nfl 
nondimens ional stress function ^cp = a^a^^ilr^ 
coefficient in expansion for i)f (see eq. (l 4 )) 

hole radius (see fig. l) 
radial coordinate (see fig. l) 

h 

an gular coordinate (see fig. l) 
nondimens ional radial coordinate (p = r/a) 

Young's modulus 
secant modulus 
shear modulus 
Poisson's ratio (elastic) 

exponent in Ramberg-Osgood equation (eq_. ( 5 )) 
stress-energy density (defined hy eq. ( 17 )) 
con^jlementary energy 
modified con^jlementary energy 

nondimens ional modified con^jlementary energy 
integers 1 , 2 , ... 10 
integers 0 , 1 , 2 , . . . 10 
integers 

Integers 0, 1, 2, . . . 
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a, 3 , 7 integers 1 , 2 

5q_o Rronecker delta (equals 1 if a = 3 , equals 0 if 

a 3) 


Superscripts: 

E elastic part; 

P plastic part 


STRESS-STRAIN RELATIONS 


On the "basis of the simple deformation theory of plasticity, the 


stress components a^., C 0 , aad Tj,q 
corresponding strain components 

= I (Or - + 

= I (‘^G - + 

^rG ~ G ' 

In equations (l), the quantity Eg 
the uniaxial stress -strain curve at 
given hy 


(see fig. l) are related to the 
6 q, and y^Q as follows: 



is defined as the secant modulus of 
an effective value of stress Cq 


cie = (2) 

The first term in each of equations (l) constitutes the elastic 
part of the strain; the second teim, in each case, is the plastic part. 
In keeping with the usiial assim^jtions of plane stress, normal and 
shearing stresses in the direction of the plate thickness are considered 
to he negligibly small. 

The sln 5 )le stress -strain relations (l) should, strictly speaking, 
he used only if Ce is continually increasing - that is, as long as no 
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unloading takes place. However, it will "be assumed in the ensuing 
calculations that equations (l) apply regardless of whether or not 0g 
decreases; the extent to which the consequent solution actually exhibits 
unloading ■vrill he examined a posteriori. 

In order to afford a complete analytical specification of the stress- 
strain law, an analytical formulation of the uniaxial stress-strain rela- 
tion (which determines Es) is desirable. One such formulation that 
appears to be useful for a variety of structural materials has been pro- 
posed by Ramberg and Osgood (ref. 2) and is given by 


e 



(3) 


where E is the elastic modiilus, is the value of stress at which 

the secant modulus Eg is equal to 0.7E, and n is a parameter chosen 
to provide the best fit to the stress-strain curve of the actual material 
^lnder consideration. Alternatively, may also be considered as an 

arbitrary parameter that may be adjusted to provide a good overall fit 
to the actual stress-strain curve, and the requirement that it specify 
the actual stress at which Eg/E = 0.7 may accordingly be dropped. 
Eqtiation (3) may be recast into the form 



(5a) 


Figure 2 shows plots of against Ee/02_ for .values of n = 5^ 5.? 

9, 19^ and 00. The gently sloping stress-strain curve for n = 5 is 
typical of some stainless steels; the sharply breaking curve for n = I9 
is similar to the' stress-strain curves of some al.umi num alloys . 

With the use of equation (5)7 the secant modulus Eg needed in the 
stress-strain relations (eqs. (l)) may be written in terms of the effec- 
tive stress 0g as 


Es 


E 
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ANALYSIS 

Elastic Solution 


Under the action of an applied stress Ooo at infinity, the elastic 
stress distribution in the plate shown in figure 1 is given by (ref. 5^ 

p. 80) 



(5) 


where p = r/a. 

7t 

At the hole (p = 1), the maxiflam value of C0 occurs at 0 = ±— 

and is equal to Thus, in the elastic range, the stress concentra- 

tion factor is 5* 


Assumed Fonn of Plastic Solution 


The plastic stress distribution may be written in the form 

CTt, = CfT'° + 


= ag- + 0Q 


’’’rQ ~ ^r0 *"^r0 


( 6 ) 
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and attention may be directed to the determination of the "correction" 
stresses Cq, and which, ■vdien added to the stresses coiresponding 

to the elastic solution, yield the stress distribiition in the plastic range. 
Since the correction stresses must be self-equilibrating, they may be 
expressed in terms of a stress function cp(r, 0) eis 



A 


1^1 



> 


(T) 


r - 1 

^r0 


^0 


.2 

1 ^ CP 
r ^r S0 

> 


By letting cp = a^a-j^\|r( p , 0 ) the correction stresses, normalized with 
respect to the nominal srleld stress d]_, may be conveniently expressed 
in terms of ilr as 


^r 1 ^ ^ 

P ^ 


''i 

= Ji. M 1 
^1 P 


(8) 


It will be assxmied that, at each value of the radial coordinate p, the 
function ijr may be expressed as a Fourier series; thus. 
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T}r = ^ f (p) cos 2g0 (9) 

q=0 


■where the Fourier coefficients f^ are functions of p. Since the 

correction stress must leave the boundary of the circular hole stress- 
free, it is necessary that ffr = = 0 at p = 1. Hence, it foUo'ws 

from equations (8) and (9) that 

^ ( 10 ) 


fq(l) = fq*(l) = 0 


(q = 1,2,5,...) (11) 


Furthermore, the correction stresses must all -vanish at infinity, and, 
indeed, they -will be assumed to vanish as l/p^. The plausibility of this 
ass ump tion stems, in part, from an examination of the elastic solution 
(see eqs. (5)), which sho-ws that the effect of the hole in perturbing 
the uniform stress field drops off as l/p^ in the elastic case. It 
seems reasonable to expect that plasticity will not Introduce stronger 
perturbation of the uniform stress state. Combining this assxnrgotion with 
equations (lO) and (U) leads, then, to the folio-wing assumed expressions 
as approximations for the f\mctions fq(p) : 



where the coefficients apq are, of course, as yet iihknown; the stress 
fmction may th\is be approximated by 
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apq^p“Pcos 2qO 




•where M and N are Integers to be chosen as large as is practicable. 
If the substitution t] = l/p is made, eq^uatlon (l4) becomes 




= / (1 - n ) ^ + (1 - ^)‘ 

^ p=0 



p=0 g=l 


Spq^TlPcos 2q_0 


(15) 


In terms of deri'vatl'ves -with respect to ti, the relationship be'tween 
stress and stress function is 


.,,3 + ^2 ^ 


^ + 2^3 ^ 

“1 Sr, } 


•re 


(16) 


,^2 ^ ^ ^5 

^9 ^0 St) 


Solution of the problem (for a given material stress- strain curve) 
no-w depends -upon the determination of the coefficients apq^ for gi'ven 

values of the applied stress at infinity. This determination is made 
possible through the use of the variational principle to be discussed 
in the next section. 


Variational Principle 

Variational principles governing the solutions to boundary-'value 
problems associated with a -variety of stress-strain la-ws of plasticity 
have been re-viewed coraprehensi-vely in reference These -wariationEil 
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principles, ho-wever, are specifically limited in applicability to problems 
involving finite domains; for this reaaon, it -will be necessary to introduce 
a new modified variational principle for the determination of stresses in 
the infinite two-dimensional region presently considered- A similar situa- 
tion in fluid dynamics was encountered by Chi -Teh Wang (refs. 5 s^^d 6) who 
found it necessary to modify Bateman's variational principle in order to 
render it applicable to con^ressible-flow problems in infinite domains. 

For deformation- type theories, the variational principles for the 
stresses Involve the so-called stress-energy density, which in the plane 
stress case can be written as 


0^0 ^'^r0 

,v(br,ff0,Tr0) = / (er dc^ + gq dog + dTr0) (l?) 


where the strains are considered as functions of the stresses; these 
functions are aasumed to be such that the line integral in equation (l7) 
is path independent- (This condition is satisfied by all, deformation 
theories that have been seriously considered, including the sin 5 )le defor- 
mation theory.) For a finite two-dimensional domain A having stresses 
prescribed over its bomdaries, the complementary energy is defined as 


F = 



w dA 


(18) 


The variational principle (ref. 4) states that, when F is written for 
the true solution to the boundary- value problem. 


&F 


= 0 


( 19 ) 


for an stress variations SCp, 6(Tg, and that satisfy eq.uilibrium 

and provide no stress resultants on the boundaries. This principle is 
readily verified by noticing that, by the use of the definition (I 7 ) 
for w. 


5F — / (ej, Sch: + eg BcTq + 7j,g 5x^.0 )dA 

^ A 


( 20 ) 
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But since the strains are coHqpatihle and the stress variations satisfy 
eq.uilihrium and provide zero boundary tractions, the principle of virtual 
work implies that the right-hand side of eq,uation (20) must vanish. 
Furthermore, if 5^w > 0 for all stress variations, 8^ > 0, and hence 
F is a relative minimum. 


y V dA = (up 6cTr + uq 8T^)ds 




Now suppose that the external boundary of A recedes to infinity in 
all directions, and that the boimdary- value problem specifies a constant 
stress state at infinity and vanishing 
tractions along all internal boundaries. 

Define Ar as the area between a large 
circle Cj^ of radius R and the inter- 
nal boundaries. (See sketch.) Then, 
by the principle of virtual work. 


\ 


/ O 


\ 

\ 




where Uj. and ug are the true displacements in the radial and circum- 
ferential directions. Hence 



““ 


5 

/ w dA - j 

^ + ugTrg)ds 

Cr 


Jar J 


( 21 ) 


for all admissible stress variations. Unfortunately, this variational 
equation is not of much practical iise as a tool for finding the stresses 
because, for any finite value of R, the displacements Uj, and ug are 

not known a priori. Furthermore, it is iiseless to let R become infinite 
inasmuch as the quantity in the brackets would then become infinite. 

The situation may be remedied by subtracting from the bracketed 
quantity of equation (21) certain integrals that are Independent of the 
plastic solution and that keep the expression finite as R becomes 
Infinite. Th\is, if the stresses of the solution to the elastic problem 
are denoted by Cr®, hg°, and Tpg°, the modified equation may be 

written as 

5 jy' [v(trr^cTe,Tre) - v(cTr°><J9°Ar9°fjdA ~ ° t °) + U9(Tr9 - dfl • 


= 0 (22) 
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Equation (22) is valid for any given value of R because the extra 
quantities in this equation, ■which depend only on the elastic solution, 
contribute nothing to the -variation. If the differences be^tween the 
stresses of the elastic and plastic solution are assumed to -vanish at 
infini-ty as l/r^ or faster, it can be sho-wn that 




lim i 
R — 


- w(cTi.°j‘f0°>Tr0O)^dA - [%(dr - ^(^re - Tr0°^ ds 

( 25 ) 




is fini-fce. Furthermore, there is no difficul-ty in e-valuating the line 
integral in equation ( 23 ) since the displacements at infinity must be 
asymptotically equi-valent -fco the easily calculable displacements that 
-wo-uld occ\ir for the infinite plate -without any interior holes. Thus, 
since equation (22) is -valid for all -values of E, it follows that 
6Finod = 0 all admissible stress -variations, where now the conditions 

of admissibili-fcy must be extended to permit only stress -variations that 
■vanish at infinity as l/r2 or faster. 

It is now possible to apply the direct methods of the -variational 
calculus (for exan^le, the Rayleigh-Eitz method) by substituting into 
equation ( 23 ) expressions for the stresses that contain undetermined 
parameters and then applying the condition that Fmod must be stationary 
■with respect to these parameters . If it is kno^wn that 6^w > 0 for all 
stress -variations, the stronger condition that F^od ^ minimum -with 

respect to the parameters can be -used; this condition for a relati-ve 
mi.ninnmi is actually satisfied by the simple deformation theory (ref. 4) . 

The general expression for Fj^q^ given by eq\aatlon (23) can be 

sinqjlified; the reduction is most economically performed -with the help 
of tensor notation and, is gl-ven in appendix A. The result is 


W= 

R-^ 


f 




^0r^ + ^0^ - o'r'^e + 5Tr02) 
2E 


dA 


- / (iir^r + U0T3.0)ds 


(24) 
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where 


w^’(ar,(Je,Tre) = ^ ^ + ee^dffg + /re^^^rej 


Uj. = Up - Up° 


ll 0 = Uq - Uq 


o 


Here Up° and Uq° denote the displacements associated with the elastic 

solution to the problem- The superscript P on the strains denotes the 
plastic part; that is, CqP, and 7^9^ are given by the second 

term in each of equations (l), respectively. 

The expression for Fmod reduced to the form appropriate 

to the present hole problem and the simple deformation theory. Prom 
equation (2) it follows that 


or 


cTe Ic^e = (°’r - 2 


- ap] da0 + 5Tp0 dTp0 


(25) 


Hence, by eqmtions (l) and (2), 


e/dcTp + €0PdaQ + rre^'i'^re = " iVe 


Es E J 


so that, with the use of equation (4) for Es, 

/ 

P P P 5 Pe \ CTe 

£p dcTp + €0 dff0 + 7p0 dTp0 = ylo^l "e" 
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Conseq-uently^ 


-w^(cJj.,aQ,T^Q) = 





whence 


wP(ffr^cr0^Tre) 


3 

7(n + 1) 




(26) 


The line integral in eq^uatlon (22) may be evaluated by noting that, 
very far from the hole, to within a rigid-bocly displacement. 


u^ 


Ea(Ooo) 


Uy ~ -y 


+ h 


E 


1 1 \^ 

2lEs(crco) " e) “ 


where x and y are rectangular coordinates in the direction of, and 
perpendicular to, CTm. 

Hence, 



and, by a change of coordinates. 
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or 




Yi + 2 


cos 20 a„ 


Ti 

U0 ~ — — - 


,n-l> 


7E \fl/ 


— sin 20 cr„ 


( 2 ?) 


With the substitution of equations (27) and (26) into equation (24), 
the modified con 5 )lementaa:*y energy, expression appropriate to the present 
problem may be -written nondin^nsionally 


2FmodE^ r” 

J 1 Jo 


2 2 
a^a]_^ 


6 K 


n+l 


n+1 


v2 


7(n + 1) y 7(n + 1) Iffj 


p dp d0 - 


_n p 2rt 
lim / p‘ 

1^ 0 


j(l + 3 cos 20) - 5 


/- \ 

' 


|sin 20 


d0 ( 28 ) 


■where ~K = g„/<^2. P ~ r/a. 

The quantities Oq° , eind Gq are gi-ven by equation (2) -with 

the use, respectively, of the actual stress components, the stress com^ 
ponents from the elastic solutions, and the correction stress components. 
The final and most con-venient form of the energy expression is obtained 
from equation ( 28 ) by using the fact that the integration need be carried 
out only over one quadrant of the plate and by introducing the substitu- 
tion Ti = l/p. Then, the function to be minimized becomes 



lim 

T] >0 



^^^(l + 3 cos 20) 



(29) 
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It is of interest to note that this expression and, hence, the solu- 
tion for the stresses do not depend on the elastic Poisson's ratio. 


Application of Rayleigh-Ritz Procedure 

Analytical minimization .- An approximate analytical solution for 
the stress distribution could, in principle, be effected by retaining a 
specified number of undetermined coefficients in the stress function (15)> 
calculating the corresponding forms of the correction stress from equa- 
tions (l6), substituting in the expression (29) for <I>, and then minimizing 
^ with respect to the coefficients. Such a process would then yield a 
finite number of simultaneous nonlinear algebraic equations for the unknown 
coefficients. In practice, the analytical evaluation of $ constitutes a 
very laborious calculation because of the (n + 1) powers of that 

appear. The work required increases very sharply with increasing values 
of n and with increases in the number of undetermined coefficients that 
are taken into account. Consequently, an analytical solvrtlon was carried 
out only for the case n = 3 — corresponding to a gently sloping stress- 
strain curve — and with only three undetermined coefficients. 

The contribution to <5 due to the term containing the single inte- 
gration with respect to 9 in equation (29) is very easily evaluated, 
and turns out to be (see appendix B) 

■^3~(f00 5 aoi) ^50) 

This result is valid regardless of how many coefficients are taken 

into account and for all values of n. The calculation of the double 
integral for n = 5 with only s-qq^ ®’ 01^ ^"10 considered in 

the expression (15) for the stress function was, on the other hand, tedious 
and time consuming. The three simultaneous nonlinear equations obtained 

do do do . 

from the conditions = = = 0 applied to the form of 0 

^^00 ^®'01 ^10 

thus found were solved (by successive approximation) for various values 
of ~K = with the results shown in table I. 

The primary use made of these resiilts was to gain insight into the 
range of values ass\mied by the coefficients, and, becaiase of the great 
labor Involved, the analytical approach was abandoned in favor of the 
numerical process to be described in the following section. 
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Numerical minimization .- The essential idea involved in the numeri- 
cal minimization of <l) is its numerical evaluation for systematically- 
varied sets of coefficients a^g^ that lead (monotonically) to lo-wer and 

lower -\ralues of <Ii. The n\mierlcal evaluation of the double integ 2 ral in 
the expression for $ -was made as follows: The rectangular domain 

0 ^ T ^ 1 and 0 ^ 0 ^ rt/2 in the i],0 plane -was divided into a 
10 X 10 grid. With the integrand denoted by I(t], 0), the double integral 
■was then approximated by two successive applications of Simpson's rule as 
follows : 



( 51 ) 

where, according to Simpson's rule, Sq — 1, Sp - S 2 — 2, S^ — 4 . . . 
Srj = it-, Sq = 2, S^ = k, and S^^q = 1. Special treatment -was needed at 

T1 = 0, since direct numerical e-VEiluation of the integrand is not feasible 

nJt/2 

there. Fortunately, the Integral / l(O,0)d0 can be e-mluated 

Jo 

analytically, and turns out to be (see appendix B) 



l(O,0)d0 


n 


3itX 

56 


f-2a^^ 

V 00 


^ 01 * 




11 J 


(32) 


This result is veilid no matter how many coefficients a are taken into 

PQ. 

account. Combining equations (50) ^ (51)^ and (52) gives the final expres- 
sion for numerical evaluation: 


rt 10 10 

IBOO ^ ^ 




n 


jtA 

560 


( 5®^00 


+ -H 


2a + 5a 
10 ^ 11 , 


(53) 


With the use of ejqjression (55) for approximate numerical minimiza- 
tion -was performed on the SEAC. The work was limited to the determination 
of only the four coefficients 8 -qq^ ®'01^ ®'ll* ^ detailed 

description of the numerical procedure — essentially an application of a 
"method of steepest descent" — is given in appendix C. As has been stated, 
the method involves the repeated calculation of $ for systematically 
varied sets of the coefficients apg that lead to a minimi mi -value of $. 

Results were found for four stress-strain curves, described by values of 
3, 5^ 9} arid 19 for the Ramberg-Osgood parameter n; for each n, several 
values of applied stress, as specified by A = were considered. 
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The resiilts foimd for the coefficients Spg are given in cable II. Con- 
vergence to the final results was generally slow; from 12 to 35 cycles of 
iteration were made for various values of A. In each case, the final 
values of the coefficients a^eed with those of several preceding itera- 
tions to within approximately 1 percent of the largest coefficient, hut 
there was no absolute guarantee that many additional iterations would 
not have changed the results appreciably. 


RESUITS AND DISCUSSION 

•V 

Stress Concentration Factor 


Once the values of the coefficients apg^ are known, the stress state 

at any point can be found from equations (5), (6), and (l6) . In terms of 
the four coefficients considered, the stresses are then given by 


^ = |[(i - + (i - + 3n^) cos 2ej - - Ti)aQQ + (n - 

(4 - lOii + cos 20 + (511 - 12 tj^ + cos 29 j 

^ =‘| (1 + 42) - (1 + cos 29j + ti2|(i - 2ii)aoo + (2n - 

(4 ti - 6ri2)a^^ cos 29 + (2 t| - 12i]2 + 12Tj^)a^ ^ cos 2^ 

=-^(1 + 2t]2 - 5ti^) sin 20 - - Stj + sin 29 + 

(4 t] - 12ti2 + 8T]^)a^ sin 2^ 


Of perhaps the greatest Interest is the value of the stress concentration 
factor determined by the ratio of 0g(l,jt/2) to the applied stress. [This 
factor is found to be 
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K = 5 - 


apo + 2aQ2 + aiQ + 2an 

A 


(35) 


Figure 3 shows the variation of the stress concentration factor with 
A = cToo/'^l each of the stress-strain curves considered. Also 

Included in the figure is a limiting cuirve for the case n = cor- 
responding to an elastic — Ideally plastic material (fig. 2). This cirrve 
was obtained sin 5 )ly by assuming that the maximum stress at the hole is 


3a„ imtil 


Coo reaches 


and thereafter remains at the value 

(or for i < A < 1), the stress concentration 

factor is K = 3A- Values of Tv greater than 1 are meaningless for 
this case Inasmuch as the applied stress can never exceed It is 


thus, for < Coo < 


interesting to note that the curves calculated on SEAC for n = 5^ 5^ 9^ 
and 19 , as well as the limiting curve for n = <», all intersect in a very 
small interval eiround Tv = 0.^1-. 


It is of particular interest to cong)are these results for stress 
concentration factors with those obtained from the following formula 
si:iggested by Stowell (ref. 1): 


K = 1 + 


Eg(bc„) 


(36) 


For the Ramberg-Osgood stress-strain curve, Stowell' s formula for K 
becomes 


K = 1 + 


2(1 ^ 
1 + 


(37) 


The variation of K with Tv = obtained from equation (37) is 

given in figure h for each of the four stress-strain cTurves considered. 
The results of the present theoretical calculations and, in addition, 
the results obtained from the three-coefficient analytical solution for 
n = 3 shown for con^jarison. 

It is seen from figure that the agreement of the present results 
with those obtained from Stowell' s formula is only fair. Perhaps a more 
meaningful con^jarison is alforded by figure 5^ which shows the variation 
with of the mRYlTmnn stress Itself (nondimenslonallzed with the 

stress for each of the four cases. It is interesting to note that 

as n increases a tendency for the curve to develop an inflection grows 
and. Indeed, this inflection becomes so distinct at n = I 9 as to indi- 
cate a very rnnal 1 reduction of stress with Increasing load, at about 
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= 0.7- Such a reduction of stress, or "unloading," really invali- 
dates the use of the stress-strain relations (l) hixt the nagnitude of the 
\mloadlng is so small (at least at the point of the plate presently con- 
sidered) as to he probably not too important. 

The results for the stress plotted according to Stowell's formula 
deviate from that of the present theory by amounts varying up to about 
15 percent. It is in^jossible, however, to ascribe much significance to 
any agreement, or the lack thereof, between the two sets of results. 

V/hlle the present results stem from an approximate solution that satisfies 
equilibrium of stress exactly and compatibi3J.ty of strain approximately, 
Stowell's formula is based on a treatment that ignores con^jatibility 
entirely and satisfies equilibrium in some average fashion over the entire 
region exterior to the hole. The present analysis clearly has a much more 
rigorous theoretical foundation, but its actual physical validity is not 
known inasmuch as it is based on an axbitraxy plasticity stress-strain 
law as well as being approximate, even within the framework of the assimied 
theory. On the other hand, in reference 1 Stowell exhibits good agree- 
ment between the prediction of his formula and experimental results for 
stress concentration factors obtained by Griffith (ref. 7) 2024-T 

(formerly 24S-T) alimiinum alloy. A careful study of the stress-strain 
curve measured by' Griffith for his specimen reveals that it cannot be 
described satisfactorily by a Ramberg-Osgood equation; consequently, no 
meaningful comparison between Griffith's results and the present analysis 
can be made. Certainly, Stowell's formula has much to recommend it by 
virtue of its simplicity; but the ultimate assessment of its validity, 
as well as that of the present analysis, must come from experiments on a 
variety of materials having both gently sloping and sharply breaking 
stress-strain curves. 


Stress Distribution 

Distributions of stress have been computed from equations (34) and 
are shown in figure 6 for the case n = 9 “ 0*9^ these dis- 

tributions are typical of those that occur for other cases in which sub- 
stantial plastic flow occurs in the nei^borhood of the hole. The sepa- 
rate sketches in figure 6 show the variations of the stress components 
along three radial lines; 9 = jt/2 (the location of maxiimim tension at the 
hole), 0 = ot/4, and 9=0. For comparison, the elastic stress distribu- 
tions are also shown. The largest deviations between the two occirr, as is 
to be expected, -vdiere the elastically computed stresses are highest. In 
addition, a rather striking difference between the elastic eind plastic 
distributions of the circumferential stress Oq occurs at 0 = jt/4; 
in the plastic case, the circximferential stress rises to a maximum and 
then drops, in contrast to the monotonic decrease in the elastic case. 
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The variation of Og, as given by equations (2) and (54), with 

Increasing applied stress at infinity has been subjected to numerical 
study in order to detennlne ■vdiether the solution Indicates imloadlng any- 
where in the sheet . The very sli^t amount of unloading that occurs at 
the hole (where Ug = 09 ) for n = I 9 and 0 = «/2 has already been 
noted in figure 5 ; spot caTculations of Og elsewhere in the sheet fail 
to indicate unloading except at values of Og well below for which 

the behavior is still essentiaJLLy elastic. 


CCNCLUDING REMARKS 


The theoretical solution presented for the stresses in the plastic 
range ajrotind a circular hole in a plate subjected to uniaxial tension is 
very fax from a final answer to the problem considered. The solution is 
based on a stress-strain relation of questionable validity; and it is only 
approximate, even within the framework of the postulated theory. However, 
it is felt that the solution has intrinsic theoretical Interest since, 
except for problems with radial symmetry, little attention has hitherto 
been directed at plane stress problems for strain-hardening materials. 
Further, the solution fvilfills to a limited extent the promise of varia- 
tional principles as a useful tool in the solution of boundary- value 
problems of plasticity. On the negative side of the ledger is the fact 
that the mmierical minimization process used in conjunction with the 
variational principle is not very efficient. Extension of the present 
approach to include many more degrees of freedom woiild sorely tax the 
capacity of even the largest high-speed computing machine presently 
available, and woiild probably require prohibitive amounts of machine 
time. The development of more efficient numerical minimization proce- 
dvires wovild be a boon to the use of the variational approach to nonllneax 
problems of the kind considered in this paper. 

Of great interest would be another treatment of the problem via 
the simple flow, or IncrementaJ., theory of plasticity, and the subsequent 
comparison with the present results. Although the flow and deformation 
theories miay differ substantially from one another for arbitrary stress 
paths, little is known about the effect of such deviations on the stress 
distributions that would be predicted by the two theories in problems 
such as the presently considered one. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 16, 1955 • 
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APPENDIX A 

SIMPLUTCATION OF MODIFIED COMPLEMENTART ENERGT EXPRESSION 


In tensor notation, with the use of the summation convention, the 
complementary energy density for plane stress may be written 





where the indices take on only the values 1 and 2. The modified com- 
plementary energy expression (25) is, in tensor form. 



where Tq, denotes the traction on Or, is the displacement vector, 

and the stperscript o denotes the elastic solution. Introducing the 
notations hap = O’op " ^ = '"a " ‘ 


Fmod= ^ ^ 

r 


/ Ta(^° + 

R >00 


°cxp 

Ce ^ 




/ 




+ v)h^ - 
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and 


. P = 3/-l _ iVa 
ap 2\Eg eA< 


_ ^ g \ 

o-P " 3 77 apy 


Now 


r'^aP '*’^aP E , l(7, \ o- 

= eL<^ + v)0^p- 0„p 

“<I^ 






SP + *')“apV - ''”aa“^p] 


But 


g (1 + „ - V0 0„„1 = 0_.„e. 


o- 

■^ap *^ap ‘'“oa 


. ‘^pp] “ 


ap ap 


and, by the principle of virtual work. 


(A2) 


f “dpVp° ^ 

•^R ^R 


ds 


Hence, 


^mod = ^ 


^R 


'®ap°«ctp 


^ap 


e „ d 0 Q + 

aP aP 


i[(i + v)5^p5^p - vS^Spp] 


dA - 



(A3) 


Now 


(l + v)0„Q0 Q - V0 
' ' ap otp oa 


'^PP " i(' 


Cf „ - — 0 0„„ 

aP aP 5 oa pp 


) + (v - i)( 


aP aP oa 


.''pp) 


The first of these two terms is precisely as defined by equa- 
tion (2), and the quantity T • ? 

-J n 


aP P 

^aP '^’^aP w , as defined immediately 
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after equation (24) . Hence, equation (A5) is the tensor equlvaJ.ent of 
equation (24), with the additional term 


1 im 


V - i 

^ P / - — 

'5 - a )dA 






)■ 


(A4) 


It remains only to show that this term vanishes to estahllsh the validity 
of equation (24). 

In terms of the stress fimction q). 


0_=cp 6„-cp 

ap ^,77 aP ^,aP 


where, as is us\ial in tensor notation, commas denote differentiation. 
Then 


U'' 


'ap°^ap “ *^oa,^ 




^Ri- 


dA 


(A5) 


Now since 0 ^^ provides zero resultant stress at all Interior boundaries, 

cp is a single-valued function (see ref. 5> P* 191) > hence. Green's 

,a 

theorem can be applied to the area integral (A5), with the result that it 
equals the line integral 




where the line F Includes both the circle Cp and the internal boimd- 
aries, and np is the exterior unit normal to F. But 
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1 

= - / <P 5 qIIq ds 
Jp ^,a op P 

= - / cp T ds 

Jp > ^ 

Since T„ vanishes at internal boimdaries, 

Cw 

% = - r ^P,a^a 

JCr 

But since cr^p, by assumption, dies oxit at infinity as l/r^ or faster, 
it follows that 


lim Lp = 0 
R — ><x> 

Thus, it is established that the expression (Ah) is zero and, therefore, 
the sin^lifled form of given by equation (24) is equivalent to 

equation (23) • 
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APPENDIX B 


EVAJUATION OF TWO LINE DffEGRAIS 


The term 


- 11m ^ f ^ ^(l + 3 cos 20)- 3fl^i]sln 29 

Jo ^ W/ 


^de 


(Bl) 


appearing In the expression (29) for 0 is 2readily evaluated. From equa- 
tions (3^) which give the stresses in terms of the coefficient Sp^ in 

the stress fimctlon, it is seen that 


i = -l^(aoo ^01 


^ = < 


2a^p sin 2e) + o(ti5) 


The qua n tity (Bl) is therefore 

kii 

0 


fr/n + cos 29)(i + 3 cos 29) - 6aQ2_(sln 29)^Jd9 = 


3 rtX 


n 


28“(^00 5aoi) 


The evaluation of the integral 
involved. From equation (29)7 



is a little more 
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Since the correction stresses are 0 (t]?), then is n-n<^ 

therefore the last term in the integrand makes no contribution in the limit 
as T) approaches zero. The first two terms in the integrand can be 


expanded into 



- 2 

Now, terms of the order will make no contribution to the integral 

as T) — > 0 because is 0 (t] 2 ) . Hence, -using the binomial expansion 

formula leads to 
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Th\is, 




6 



de 


Now cr □ is 0(l); f lorthermore , terms in the stress function containing 

"0^*2 

Spq^ lead to terms in the correction stresses of the order . Hence, 


consideration need he limited only to p = 0 and 1. Furthermore, 


■^aP 


has only terms that axe either independent of 0, or contain cos 20 or 
sin 20. Hence, by the orthogonality properties of these trigonometric 
functions in the range (0,it/2), only values of q_ = 0 and 1 provide non- 
vanishing results for the integrals of products of the elastic stress and 
the correction stresses. Thus, no matter how many coefficients a^_ are 


taken into account, only ^qq, ®-10^ ^•ll possibly contrib- 
ute to the integraJ. being evaluated. Using the appropriate parts of equa- 
tions ( 54 ) for and and legitimately ignoring terms of the 

order or hi^er in their products, yields 


prt /2 p «/2 

/ l(O,0)d0 = lim / ±— 

^0 Ti^O^O ' 


'■00 


(- t) ®io(a-) 


aoi(6ii^ cos^ 20 - 6 t]^ sln^ 20 + 21 t]^ cos^ 20 - 

24 ti^ sin^ 20) + a^^- cos^ 20 + 12t]^ sin^ 20) 


-4- d0 


Integrating with respect to 0 and then letting t) — ^ 0 gives 


n 


yCK 

56 


(-^00 - + ^10 + 5^11) 
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APPENDIX C 

DESCRIPTION OF NUMERICAL MINIMIZATION PROCEDURE 


The numerlcaJ. minimization procedure used was essentially the same 
as that used in reference 8 in the solution of a different plasticity 
problem; much of the ensuing description parallels that contained in 
appendix C of reference 8, but is included herein for the sake of com- 
pleteness of the present report. 


The so-called steepest -descent procediire (see, for example, refs. 9 
and 10) which formed the basis of the numerical minimization may be 
described in general terms as follows: Consider a function il>^x^,X 2 , . . .x^^. 

The set of n Independent parameters may be conveniently denoted by the 


n-con 5 )onent vector x. (i = 1, 2, 3) 


n). 


minimizes ^ is sought 
the gradient of 4, that is. 


An initial trial vector 

8$ 


The value of 

. ( 0 ) 




that 


tlon - 


8f 


8xj^ 


is calculated at x^ 


is assumed, and 
(0). The direc- 


8xj^ 

the function 


is then the direction of steepest descent of the function 4>; 


<5 


values of 

( 0 ) 


( 0 ) 


- 5 


So 


8x^ 




is then evalviated for various positive 
5 that minimizes 
When this value of 8 is found (presijmably, approxl- 


5 in an effort to find the value 
8x. 


mately), a new direction of steepest descent is determined by evaluation 


of the greuilent of at the point x. 


(1) _ ^{0) 

- ^ 


5^ 

8xj^ 


I ( 0)1 


The 


process is contln\ied until satisfactory convergence is obtained to the 
lowest possible value of $ . 


In the present problem, the function $ is given by equation (35), 
and the four coefficients ®'01^ ®'10^ ®’ll role of 

the components of the vector x^. The basic procedure outlined was mod- 
ified in several respects. The evaluation of the greidient of <Ii was 
actually performed on the basis of a finite-difference approximation to 
each of the four partial derivatives required. Also, the value of 8 in 
any given cycle was found as the minimum point of a parabola through the 

where 8 was 

8xiy 

chosen to be of some convenient magnitude, preferably of the order expected 


three points 5>(Xij, <Ii fx^ 


<t 






Xi - 2 
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for 6. The procediire was fiirther modified by Inserting after each two 
successive cycles of minimization. In the direction of the negative 
gradient, a third cycle of a different nature. In this extra cycle, 
minimization was performed in the direction determined hy the difference 
between the last-obtained approximations to the unknowns and the approxi- 
mations of two cycles before. The motivation for this exbra cycle stems 
from the fact that the nnmodified method of steepest descent often tends 
to furnish successive approximations to the minimizing vector that zigzag 
toward the true minimum; the extra cycle was an attempt to speed up con- 
vergence by moving in the direction determined by the mean of a zig and 

a Z8Lg. 

The procedure outlined was coded for calculation on SEAC; once values 
of n, A, and initial estimates for the coefficients a^g^ were put into 

the machine, the iteration process proceeded automatically. Successive 
approximations to the minimizing coefficients were printed; the number of 
cycles of iteration performed varied from 12 to 35* In each case, the 
final resiilts agreed with those of several preceding iterations to within 
approximately 1 percent of the largest coefficient. Although the natiire 
of the mnnerical procedure is such that the accuracy of the final results 
can not be positively assessed, some confidence in their accuracy is lent 
by the smooth variations with A obtained for the stress concentration 
factor . 
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TABLE I 

COEFFICmraS ANALYTICAL SOLUTION, n = 5 


A 

^00 

^01 

^10 

0.2 

-0.00011 

o.ooit58 

0.00991 

,5 

.00521 

.04647 

.09035 

.7 

.01868 

.09274 

.16796 

1.0 

.05330 

.175^6 

.28840 
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TABIE II 

COEFFICIENTS NUMERICAL SOEOTION 


"00 


"01 


^10 


"11 


n = 5 


0.0103 

0.0146 

-0.0027 

-0.0109 


.0079 

.0152 

.0455 

.0244 


. 0304 - 

.0688 

.1457 

.0445 

1 

.0666 

-1551 

.2666 

.0338 

n = 5 


0.0013 

0.0004 

0.0014 

0.0017 


.0015 

.0052 

.0542 

.0436 


.0703 

.1221 

.1986 

.0652 

1 

.1572 

.2808 

.2885 

-.0250 


.2287 

.4489 

.4013 

-.1364 

n = 9 


- 0.0014 

-0.0018 

0.0082 

0.0102 


-.0154 

-.0006 

.1437 

.1159 


.1131 

.1712 

.2401 

.0764 

1 

.3549 

.5130 

.2548 

-.1866 


.4928 

.8190 

.2871 

-.4256 

n = 19 


-0.0116 

-0.0158 

0.0545 

0.0672 


-.0286 

-.0271 

.1556 

.1647 


.1617 

.1895 

.2610 

.0995 

1 

.6440 

1.1120 

-.0240 

-.6039 
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Figure 2.- Nondimens lonal Ramberg-Osgood stress-strain curves. 
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Figure 3*- Variation of stress concentration factor with 







